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In this paper the tachyonic resonance preheating generated from the bosonic trilinear (f>x 2 inter- 
actions in an expanding Universe is studied. In A<^> 4 /4 inflationary model the trilinear interaction, 
in contrast to the four- legs (f> 2 X 2 j breaks the conformal symmetry explicitly and the resonant source 
term becomes non-periodic, making the Floquet theorem inapplicable. We find that the occupation 
number of the produced x _ P ar ticles has a non-linear exponential growth with exponent ~ x 3 ^ 2 , 
where x is the conformal time. This should be contrasted with preheating from a periodic resonant 
source, arising for example from the four-legs 4> 2 X 2 interaction, where the occupation number has a 
] linear exponential growth. We present an analytic method to compute the interference term com- 

f*^ , ing from phases accumulated in non-tachyonic scattering regions and show that the effects of the 

f"-; ■ interference term causes ripples on x 3 ^ 2 curve, a result which is confirmed by numerical analysis. 

| Studying the effects of back-reaction of the x-particles, we show that tachyonic resonance preheating 

in our model can last long enough to transfer most of the energy from the background inflation field 
4>, providing an efficient model for preheating in the chaotic inflation models. 
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INTRODUCTION 



At the end of inflation the Universe is extremely cold and all of its energy is concentrated in the inflaton field <f>. 
' Therefore, the reheating stage is a key ingredient to connect a successful inflationary stage into hot big bang cosmology. 
In simplest models of inflation, such as chaotic inflation models, inflaton field oscillates coherently around the minimum 
of the potential with amplitude close to the Planck mass, Mp. 1 To transfer the energy from the background inflaton 
. field and reheat the Universe one has to consider its coupling to other fields, such as the Standard Model fields. Due 
1 to coherent oscillations of the background inflaton field around the potential minimum, explosive particle creation can 
^\ \ happen via parametric resonance which can drag most of the energy from the inflaton field and reheat the Universe 
[J S B 0j- Preheating, the stage of explosive particle creation via parametric resonance, is a non-perturbative 
effect and the produced particles are highly non-thermal. Subsequent turbulent interactions of different modes and 
• , re-scatterings would bring an end to preheating and thermalization starts @, H, 0, H H where the p erturbative 
reheating mechanism may start [HLClI. For 

a review of preheating see [TH and references therein. 

There have been many studies of preheating via four-legs interaction g 2 <fi 2 x 2 between the inflaton field and the 
resonance field \ m expanding or flat backgrounds. One disadvantage of g 2 (j> 2 x 2 preheating channel is that the decay 
of inflaton in an expanding background is not complete and the inhomogeneous inflaton particles will dominate with 
the matter equation of state, an unsatisfactory state to end (p)reheating with. 

In order to enhance the preheating mechanism other interactions are necessary, specially once the inflaton field 
amplitude is damped while oscillating around its minimum. The simplest and the most natural of these interactions 
is the trilinear interaction a<f>x 2 where a is a constant of dimension of mass. The implications of trilinear interaction 
for the chaotic inflationary potential m 2 (f) 2 /2 was studied in [l4l |. One interesting aspect of preheating via trilinear 
interaction is that during half period of its oscillation, the inflaton field changes sign and the time-dependent frequency 
squared for the x particle creation becomes negative for small enough mode wavelengths. This tachyonic instability 



enhances the preheating mechanism significantly and in [14| was dubbed as "tachyonic resonance". The analysis of 
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|14j were mainly devoted to a flat background and the effects of tachyonic resonance in an expanding background 
was briefly considered, incorporating them "adiabatically" . The tachyonic effect of a trilinear interaction was used in 
fl5j |. m the context of cosmological moduli problem and the idea of tachyonic parametric resonance with a negative 
coupling constant was also exploited in (l6| . 

Here we consider the tachyonic resonance preheating in X(j> /4 and m 2 </> 2 /2 inflationar y p otentials in an expanding 
Universe. Our main interest will be in the A0 4 /4 theory but we extend the results of [f4( for m 2 </> 2 /2 inflationary 
model to an expanding background. 

The preheating with four-leg interaction within the X(j> /4 inflation theory was studied extensively in [17j . Due 
to the conformal invariance of the inflationary potential, the effects of expanding background can be absorbed by a 
conformal transformation and one is basically dealing with preheating in a Minkowski background. However, with the 
addition of c0x 2 interaction, the conformal invariance is broken explicitly and the methods employed in [17| should 
now be modified. Breaking of conformal invariance results in an explicit breaking of periodicity in %-modes equation: 
the amplitude of the source term in x held particle creation grows linearly with the dimcnsionless conformal time x. 
Consequently, the standard Floquet theorem does not apply. This leads to the interesting result that the number 

3/2 

density of x particles produced at momentum k, n^, grows as e cx for some constant c. This should be compared 
with the standard result of the preheating due to a periodic source term, where has a linear exponential growth 
as dictated by the Floquet theorem. 

Another interesting result of our analysis is that we present an analytic method to compute the interference term 
which accounts for the phase difference between the x fields during each period of the field oscillation. The phase 
accumulated during certain number of oscillations, which we compute, makes the x particle creation to be destructive 
for certain modes, the k 2 of which we specify. This causes some ripples on the x 3 / 2 curve in the lnnk diagram. This 
result is to be compared with the stochastic nature of the phases in the four-leg preheating scenarios [I], [l7| ■ 

While our main interest is in A</> 4 /4 inflationary theory, we also study the effect of expanding background on m 2 (f> 2 /2 
theory. As one may expect, the effect of expanding background will generally suppress the preheating efficiency. As 
we demonstrate, however, it actually enhances the preheating for some certain modes. 

As the number of x particles grow their back-reaction on the dynamics of the particle production and the (/) field 
dynamics become important, slowing down the preheating and eventually terminating it. For a successful preheating 
model it is necessary that this back reaction does not become large too early, before the energy of the inflaton field 
is completely transferred into the <j) or x particles. It is also important that in the end of preheating we remain with 
a relativistic ensemble of these particle. As in the case of m 2 (f> 2 , which has been analyzed in [14| . our analysis shows 
that the trilinear interactions seems to lead to a more efficient preheating than the four-leg case. 

The paper is organized as follows. In section [TI] we study in details the tachyonic resonance for A^ 4 /4 theory. We 
demonstrate that our analytical results agree very well with the exact numerical results. In section HTT1 we repeat the 
analysis of tachyonic resonance for m 2 (f> 2 /2 theory in an expanding background. In section [IV] we study the effects 
of x-particles back-reactions on the preheating and estimate the time preheating ends. A summary of the results are 
given in section [V] Some technical aspects of the analysis are relegated into Appendices [AUDI 



II. THE TACHYONIC RESONANCE IN A0 4 /4 THEORY 

As explained above, we are mainly interested in tachyonic resonance in \<t> /4 inflationary theory with the trilinear 
interaction atpx 2 where a is a parameter with dimension of mass. As in [I4[ we have to include the self-interaction 
A'% 4 to uplift the potential and keep it bounded from below. The total potential is 



v- i ? + - 1 tf + - i s + (1) 



A , A a 2 A' 

+ ^x + T 

a (t-4AaO + 4 x + i) • (2) 



The potential has a global minimum at (f>o = —y / a 2 /2XX' and Xo — ~ c0o/A'. The last term in is added to lift 
the global minimum to zero. In order to make sure that this constant value does not contribute to the inflationary 
dynamics, we require that er 4 /AA' 2 <C Xtfi 4 nd where 4> en d is the value of the inflaton field at the end of inflation when 
the onset of preheating starts with 4> en d — Mp/ ' \pK. This leads to 

« I . (3) 



If we assume that the %-field is heavy during inflation and is settled down to its minimum at xo this requires that 
1 2crMj,/A</) 3 > I, which can be satisfied if YlaMpj\^\ > 1 where <fii is the initial value of the inflaton field at the start 
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of inflation. For inflation to solve the flatness and the horizon problem, we require that 4>i — WMp and A ~ 10 14 . 
So the assumption that the field x would be stabilized to its minima requires that 

> 10 2 . (4) 



AM, 



However, it is also possible that the field x i s light during inflation which can then contribute to iso-curvature 
perturbations. 

To study the preheating with trilinear interaction, we review the background of [l7T ] where the parametric resonance 
with four-legs interaction for A</> 4 /4 theory was investigated. As in [17| . performing the conformal transformation 
if = a(rj)(f) the background equation for the tp field to a very good approximation is simplified to 

ip" + Xip 3 = , (5) 

where the derivatives are with respect to the conformal time drj — dt/a(t) where a(t) is the scale factor. The solution 
to this equation is given in terms of Jacobi elliptic cosine function 

(p = <pf(x) , (6) 

where (p is the amplitude of the oscillations at the onset of preheating which we take to be O.IMp, and the dimensionless 
conformal time x is defined by 

x = V\<pr) = f - p j Vt. (7) 



The Jacobi elliptic cosine function satisfies the relation f' 2 — i(l — / 4 ) with the solution denoted by 
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f{x) = cn(x, ~) . (8) 

The function f(x) is a periodic function with the periodicity 

T = AK{\/2) ~ 7.416 (9) 

where K(m) represents the complete elliptic integral of the first kind. For a review of Jacobi elliptic cosine function 
see Appendix [Aj 

One also notes that the following relations also hold 



2tt\ (p 1 rhr (p In A <p 2 

YWj 1 ~ V TMp x ' t_ VTm„ 



a(v) = \ j., v = \ ittt x ' t = \ -r-rrv • ( 10 ) 



Ignoring the back-reaction of the produced x particles, performing the conformal transformation x = a ( 7 ?)x J an d 
using the background equations (|10[) . the equation for the mode function Xk m the momentum space is 

Xk + (* 2 +pxf(x))xk = 0, (11) 

where the derivatives are now with respect to the dimensionless conformal time x given in ([7]), and 

2 k 2 [2tt ( a 



K =V' ps VTUd- (12) 

In writing (fTTj) we have ignored the xo and <fio terms, associated with the minima of the potential |T]). The effects of 
non-zero values for xo and 0o will appear in pip as the shift in p by an amount of the order p&§- (and similarly for 

the k 2 term), causing an error of the order py^, which using ([3]) it is much smaller than one. 

Before we start the analytical theory of tachyonic resonance, we need to estimate the magnitude of parameter p. 
This in turn is determined by the conditions whether the field x is heavy or light during inflation. If we assume that 



2 Our definition of Jacobi elliptic cosine conforms to that of [19H which is slightly different than that of |17|| . 
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FIG. 1: A plot of pxf(x) given by (|13[> is presented here for p = 0.05. The horizontal line represents the solutions of 
u> 2 (x) = k 2 +pxf(x) = for k 2 — 1. The filled regions shows the tachyonic regions where uj 2 (x) < 0. In our notation, the j th 
tachyonic region is confined to x J < x < ad". 

the field \ is heavy during inflation, then the condition (J4J) indicates that p > 100. However, if we assume that the 
field x is light during inflation, which is our preferred choice, then ([3|) can be used to estimate the bound on p. Writing 
([3]) in the form pWX/X' <C 1, we see that for natural choice of A ~ A', p -C 1. On the other hand, for A' 2> A, the 
bound on p may become relaxed. We present our analytical results for arbitrary value of p. However, for numerical 
examples we shall consider the cases p <C 1 and p ~ 1 for illustrations. 



Equation (|lip represents a harmonic oscillator with the time dependent frequency 

lu 2 (x) = k 2 +pxf(x) — k 2 + px cn(x, — ) . 



(13) 



A plot oi px f(x) is given in Fig. [TJ We see that lo 2 {x) is not periodic and its maximum value is increasing linearly 
with the conformal time x. Consequently, the Floquet theorem for particle creation via parametric resonance with 
a periodic source does not apply. We need to employ the direct scattering method to evaluate the number density 
of the x particle, n^. For a given k 2 ^ 0, uj 2 (x) is initially positive. However, for t > t* it becomes tachyonic and 
the method of tachyonic resonance considered in jlij would apply. Interestingly, the zero momentum case, n = 0, is 
tachyonic for all the time. We therefore start our analytical study of tachyonic resonance with this simple case and 
defer the case of k 2 ^ to the section Hi Bl 



A. Tachyonic Resonance with k = 



As can be seen from Fig. [TJ oj 2 (x) becomes negative in the region where f(x) < 0. The Jacobi elliptic cosine is 
periodic with period T = 4K(l/2) and has roots at 



3_ , ,. 1. 



*7 = C?- 5 )t • 4 = C?-i) T > ( 14 ) 

where the integer number j = 1, 2, ... counts the number of oscillations of f(x) starting from x — 0. In this convention, 
the j tachyonic region is confined to xj < x < x^ . Similarly, the minima and the maxima of f(x) are given by 



x? ax = (i - 1)T , xf n =(j--)T. (15) 

In the region xj < x < x^ , lo 2 {x) < and the method of tachyonic resonance developed in [l4| applies. 

As in [l4j, in the region uj 2 {x) > and when the adiabadicity condition /oj 2 <C 1 holds, the solution of (fTTj) can 
be given in WKB approximation as 



Xk=o(x) - 7= . . exp f-i [ cj(x')dx') + — exp ( +i [ u{x')dx' ) , 

y/2u)(x) V Jx J V 2w ( X ) V Jx o 



(16) 
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where a? and (3 J are the Bogoliubov coefficients with normalization |a| 2 — |/3| 2 = 1. At the beginning of preheating 
when the inflaton field starts its oscillation toward the potential minimum, there is no \ particle and we start with 
the vacuum initial condition a — 1 and j3° = 0. The occupation number of the x particle after j oscillation (or after 
j tachyonic regions) is given by 

4=o = l^' I 2 - (17) 

For the tachyonic region xj < x < atj, the WKB approximation holds again for the frequency f2 2 (x) = —lo 2 (x) and 
the solution is given as a superposition of exponentially growing and decaying parts: 

Xk=oO*0 - -4= exp ( - / Q(x')dx' ) + Z— exp ( + f Sl(x')dx' ) , (18) 
^2Q.(x) \ J x i_ J ^2fl(x) \ J x '_ J 

where a? and b 3 are constants of integration. Finally, after x > x 3 + , the solution is given by (fT6|) with j — > j + 1. 

Around the points xi_ and where the adiabatic approximation is broken, we have to solve the mode equation 
(fTTj) as in scattering theory and match it with the solutions (fTB| and (|18p . Following the methods of [l4[ in performing 
this matching condition, we obtain the following transfer matrix, relating the Bogoliubov coefficients after j th + 1 
scattering to those of j th scattering via 



3+1 




where 



(19) 



X 1 = I + Q(z)dz, (20) 



and 9 J is the total phase accumulated from Xq to Xj during the non-tachyonic intervals where u> 2 > 0. In our case, 
we have J = 8° + J^j where 6° is an initial phase and 

r- 

9 J = / cu(z)dz. (21) 

The key difference in our case compared to the case of m 2 (f) 2 theory in flat background studied in [l4[ is that, because 
of the non-periodicity of u 2 (x), The local maxima and minima of uj 2 (x) increases linearly with j. This results in a 
non-trivial j-dependence in X J and Q 3 . Consequently, the occupation number of the \ particles after j oscillations 
of the inflaton is 

4=o = l^l 2 = exp (£>XM f[(2cose s ) . (22) 

\l=l ) 8=1 

In order to provide some useful analytical expression for n£. =0 , we need to perform some reasonable approximations 
in calculating the integral for X 3 and <d 3 in (|2"01) and (|2~lj) . To calculate the integral in X 3 we note that the bounds 
of integration in ([20)1 ranges from a;™ m — T/4 to xj 1 " 1 + T/4. This suggests that as a good approximation, we can 
approximate O(x) = pxf(x) ~ px m i n f(x) and the integral in (f2"0"l) is approximately (for the details see Appendix [C]) 



'""+T/4 r + % 



l -T/4 



r(|) 



= / n(x)dx~^/pxf n j dx = vW/2) ^x/^T , (23) 



where r(x) is the Gamma function, and is the complete elliptic integral of the first kind [1||. Similarly, for Q 3 

c +T/4 



Q j = u{x)dx~ Jpxf ax y/f(x)dx = ^ P nK(l/2) tAjtV^J ~ 2 . (24) 



-t/4 v - j_t r^, 

Using the approximations proposed for harmonic numbers summation in Appendix |B] one finds ^Zi=i V2« — 1 
2 v / 2j 3/2 /3 and 



£ 2^ ~ | ^2pnK(l/2) j 3 / 2 . (25) 

i ^ 1 IsJ 
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FIG. 2: Occupation number lnrik=o as a function of conformal time x. In the left figure with p — 0.03, the wavy solid curve 
(blue) shows the exact numerical solution, the smooth growing curve (green) shows the analytic solution (|26[) which does not 
contain the interference term, the dash-dotted curve (red) shows the analytic solution (|27|l which includes the interference term 
added and the bottom solid curve (black) shows the behavior of the interference term. In the right figure, lnrik=o is plotted 
for different values of p. Again the solid curves (blue) show the exact numerical solutions of while the dashed curves (red) 
shows the analytic solution (|27[) . The graphs from bottom to top, respectively, correspond to p = 0.03, p — 0.5 and p — 1. The 
agreement between the full numerical results and our analytical formula, (I27|l . is impressive. 



After j oscillations, we have x — jT, so for large j the exponent for 7i k =o in (f22j) behaves as 

n k=0 oc exp(0.490v^a; 3/2 ) . (26) 

Quite interestingly, the occupation number grows exponentially with the exponent proportional to a; 3 / 2 . One can 
trace the non-linear exponential enhancement of the occupation number to the explicit breaking of the periodicity 
of u 2 (x) in (|13[) . This is in contrast to conventional model of preheating where the source term is periodic and the 
Floquet theorem applies with the result that the occupation number has a linear exponential growth. One can check 
that the breaking of periodicity in our case is a direct consequence of conformal invariance breaking via trilinear 
interaction in an expanding background. 

A plot of n k =o is shown in the left graph of Fig. [2] containing the full numerical solutions and our analytical 
results. Although the non-linear exponential growth captures very well the overall behavior of n k= o, but the figure 
from the exact numerical solution shows some small wiggles and ripples which can not be captured by the exponential 
profile. To take these minor but interesting discrepancies into account, we should also add the effects of the phase 
term originated from the products of cosines in (|22[) . As one can see in Fig. [2] this term is oscillatory but not 
periodic which can explain the small oscillations in the profile of lnnk=o versus x. Taking the effect of the phase term 
into account, one obtains 



<=o - ex P 



j 

0.490yp x 3/2 In (4 cos 2 (7.427V5^)) 



(27) 



where j — x/T. 

In the right graph of Fig. [2] we have plotted n k= o for different values of p. As this graph shows, with the addition 
of the effects of the phases, our analytical result l|27p shows a perfect agreement with the exact numerical results, 
confirming that our analytical solution captures the correct p-dependence. As one expects, the larger the value of p 
is, the stronger is tachyonic resonance from trilinear interaction. 



B. Tachyonic Resonance with k 2 ^ 

After presenting the analysis for the simple zero-momentum case, in this section we present the analysis for arbitrary 
momentum, k 2 7^ 0. As can be seen from Fig. [T] for k 2 7^ there is no tachyonic resonance from the beginning 
where n 2 > px. After some oscillations co 2 (x) becomes negative for x>p/n 2 and our previous methods for tachyonic 
resonance would apply. From (fT3"|) , we see that the frequency of oscillations becomes tachyonic after j — oscillations 
where 



(28) 
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Here [z] represents the integer part of z. 

For small x there is no tachyonic region from the start of preheating, but if the frequency becomes non-adiabatic one 
should expect particle production via parametric resonance as in conventional preheating analysis. However, as we 
shall see from our full numerical results, the particle creation from the non-tachyonic scattering for j < oscillations 
is quite negligible compared to tachyonic resonance particle creation after j > oscillations. The non-adiabadicity 
condition \lj'(x)\/u)(x) 2 ^ 1 is satisfied for 



> 



0.8472 



Jnad — 



ITy/pT 



J* - J* 



OA 

VP 



(29) 



This indicates that j na d — j* ~ few and the onset of parametric resonance is actually around j* where tachyonic 
resonance starts. Besides the very short period of parametric resonance, the "effective" Floquet index during this 
period is much smaller than that of the tachyonic regions. Therefore, the particle production is mainly concentrated 
in the tachyonic regions. This result is also verified in our numerical investigations. 

Repeating the same methods as in the previous section, for j > j*, the Bogoliubov coefficients before and after j th 
oscillations are related by 



< +1 
k +1 



,X2 




(30) 



where 



px'- 



f(z)\dz, 



(31) 



with rf in 



k 2 jpx™ ,m and Qi is the total phase accumulated from xq to x ■ during the intervals where uj 2 > 0. We 



have 6 J = 6»° + E, © 3 wi th 



QJ 



px 1 - 



+ 



f(z)dz, 



(32) 



where 



„max — ,.2 



K 2 /px 



rna'j 



Similar to the analysis resulting in (|2"5|) and using the formulas presented in Appendix [D] one obtains 



E 2 ^ 

j.+i 



n(2coseo 2 

i.+l 



1 - 2b' k 2 



1 



pTy/W^l 



4a 

y 

= exp 



^ In 4 cos 2 ay^pTy/j - 1 
j'.+i 



bK 2 



/pTVJ^T 



(33) 
(34) 



where (|B2I) and (|B3[) have been used and 

a = 2.72, b' = 2.86, 6 = 3.75 
After j th oscillations, x = jT, and the occupation number is given by 

Wk 2 



exp 



f^Oy/P t 
V 3T V 



3/2 3/2 
X ' — xJ 



\/pT 



12 - x 1/2 



x Interference term , 



(35) 



(36) 



where x„ = jJT and the interference term is the expression (|34| . 

As in the zero-momentum example, the occupation number has a non-linear exponential dependence with the 
leading exponents being a; 3 / 2 and x 1 / 2 respectively. The fact that before the tachyonic regions, i.e. for x < a;*, there 
is no particle creation is clearly seen in (|36|) . Both of these two non-trivial results are a consequence of the violation 
of the periodicity in the resonance source term. The left plot in Fig. [3] compares our analytical result (|36|) with the 
full numerical results. The agreement is again impressive. 
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FIG. 3: The left figure shows lnrik as a function of x for p = 0.05 and k 2 = 1. The wavy solid (blue) curve shows the exact 
numerical solution of {TTJ. The (red) dashed curve shows the analytical solution (|36[l . The smooth (green) solid curve shows 
the analytical solution without taking into account the interference term. The bottom (black) solid curve shows the effect of 
the interference term (jMJ. The agreement between the full numerical result and the analytical formula (|36jl is impressive. In 
the right figure shows the destructive effects of the interference term in (136[1 . The lower and upper wavy (blue) solid curves 
show the exact numerical results for k 2 = 2 and n 2 — 1.5, respectively, with p = 0.05. The lower and upper dashed (red) curves 
show the effects of the interference term, for k 2 — 2 and k 2 = 1.5, respectively. We see that the interference term causes a 
suppression of particle production for k 2 = 1.5. 



The interference term in (|36[) has a very interesting effect. Suppose for the moment that we do not take into account 
the interference term in (|36[) . Since b 1 > 0, from the exponential term one may expect that the larger the value of k 2 , 
the more suppressed will the particle creation be. Our numerical investigations supports this general rule, but there 
are some noticeable exceptions. In the right graph of Fig. [3] we have presented an example where a larger k 2 has a 
higher occupation number than a smaller k 2 . The resolution to this apparent paradox relies on the destructive effects 
of the the interference term. Looking into the form of the interference term (|34[) one observes that the interference 
becomes destructive in a region in k 2 — j space where cos(0£) vanishes. For the latter, one requires the phase term 
0£ to be stationary on its roots, i.e. cos(0£) = djO^. — at some stationary time j s . From (f3"4"| one easily finds that 
the stationary points of occur at 

bn 2 , . 

.h = 1 + -= , (37) 
apl 



with 



9£ = 2Vabn . (38) 



As explained above, demanding that 0^ — > (2n — 1)§ for some integer n, one obtains the value of K where the 
destructive feature of the interference term is pronounced 

(2n-l)7T , 

K sta b = - ~ 0.25(2n- 1) , n = l,2.... (39) 



Wab 

Note that 0^ s and K s tab are independent of p. Note also that the stable bands are equally-spaced in the mo- 
mentum space. It is worth checking this equation for some n. For example for n = 2, 3, 4 and 5, one finds 
k 2 = 0.54, 1.51, 2.96and4.89 respectively. Figj4] shows behavior of 0£ as a function of number of oscillations for 
three given values of k 2 . As is seen in the left and right figures the stationary point of can be near (2n — l)ir/2 
which causes particle production suppression. 

One also observes that the destructive phase in the interference term is not static in time. This is somewhat similar 
to stochastic resonance phenomena observed in resonance preheating in an expanding background [l|. However, 
in contrast to the stochastic resonance, in our case the interference term is not completely stochastic and one can 
keep track of the phase at each oscillation period. The value of k around K s tab, where the destructive effects of the 
interference term is significant, form semi stability bands. These semi stability bands, as indicated in (|37p . occur at 

7T 2 (2n— l) 2 

^stab = T + (2n - l) 2 ~ 7.416 + 0.0834^ . (40) 

loa z p p 
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FIG. 4: Phase accumulated through j oscillations for different values of re 2 , left: re 2 = 1.51, middle: re 2 = 2 and right: re 2 = 2.96. 
The angular position of each (red) dot indicates the value of Q J k . The outer points correspond to smaller value of j and as j 
increases the dots become denser as can be seen from the form of 0£ given in (134[) . As Gj, tends to (2n— 1)tt/2, the contribution 
of j th oscillation in the interference-term suppresses the particle production. For re = 1.51 (left) the stationary point of Qi is 
near 7r/2 . For re = 2 (middle) the stationary point is far away from the poles and for re = 2.96 (right) the stationary point is 
near 3n/2. This explains why re 2 = 1.51 and re 2 = 2.96 belong to the "semi-stability bands". Note that ln|2cos0^| is positive 
for |0{ — nn\ < n/3 and is negative for \Q J k — (2n + 1)tt/2| < 7r/6, which are separated in this figure by + and — signs. 



100 




j X 

FIG. 5: In the left figure the contour plot of the phase term in re 2 — j plane is shown. The darker regions show the semi-stability 
bands in which the interference term suppresses particle creation. The right bar represents In |2 cos Q 3 1 as a function of number 
of oscillations j. An important point of our analysis is that the semi-stability regions are independent of p. In the right figure, 
In rik as a function of conformal time for the predicted semi-stability band re 2 ~ 3 is shown. The (blue) solid line shows the full 
numerical solutions and the (red) dashed curve shows our analytical result (I36|l . As we predicted, particle production for this 
momentum is suppressed for small times because of the interference term. 

A plot of the effects of the interference terms is shown in Figj5]with p = 0.05. As an example of the evolution of 
the semi-stability band consider the line n 2 = 1.5 which corresponds to n = 3 in (f3"9"| . For small x, the interference 
term has a destructive contribution. As the time goes by, effects of the interference term would become milder and 
less important. As explained above, for this reason we may call k 2 = 1.5 a semi-stable band rather than a stable 
band. As shown above, the other semi-stable band is at re 2 ~ 2.96. The tachyonic resonance for this case would 
start at x* — j*T ~ 60. Nonetheless, as depicted in the right graph in Figj5] particle creation is negligible until 
x ~ 110; compatible with a; s t a b — 104 obtained from (|4T)|) for n = 4. This is a consequence of the semi-stability band 
for re 2 = 2.96 which suppresses the particle creation at early time and eventually loses its importance as the Universe 
expands further. 

C. Expanding vs. non-expanding background 

In the last two subsections the effects of expanding background were included in the preheating analysis via the 
conformal transformations ip — a(rj)4> and x = a ( r ?)x- One consequence of the expanding background was the 
appearance of the non-periodic factor x in lu(x). This in turn leads to a non- linear exponential particle creation. This 
should be compared with preheating with four legs interaction, g 2 <f> 2 x 2 , in Xtfi 4 /4 theory which preserves the conformal 
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X 

FIG. 6: Comparison between particle creations in a flat and expanding backgrounds in A</> 4 /4 theory. For p = 0.1 and k 2 = 0.54 
which is in tachyonic region in a flat background, but because of effect of interference term (|43[) there is no particle production 
for these parameters. Nonetheless, due to non-static interference-term, there is particle production in expanding universe. 
The wavy (blue) solid curve indicates the full numerical solution for an expanding Universe, the (red) dashed curve shows our 
analytical solution and the straight (green) line at the bottom represents the particle creation in a flat background. 



invariance and uj 2 (x) is periodic and one obtains a linear exponential particle creation. One may ask what would be 
the situation if one considers preheating with trilinear interaction in a non-expanding flat background. The trilinear 
term in a flat background is periodic and, as a consequence of the Floquet theorem, one obtains a linear exponential 
particle creation. This seems paradoxical, noting that the expansion of the Universe usually suppresses the particle 
creation. As we show below, the root of this apparent paradox relies on difference between conformal time x and the 
cosmic comoving time t. 

In a flat background with a(t) — 1, the solution to the inflaton field is 4> — 4>f(z) where 4> is the initial amplitude 
of the inflaton field at the start of preheating and f(z) is the Jacobi elliptic cosine function. Here we defined the 
dimensionless times z = yf\<t>t. The equation of motion for the resonant field x is 



^ + ( K 2 +px f( z )) X]i = Q (41) 



with x = y/3/2Tr(Mp/ip) and k 2 = k 2 /\(p 2 . 

Following the same methods as in the last two subsections, the occupation number is 

n{ = exp(2jX fc )(2cose fc ) 2 ^- 1 ) (42) 

with the crucial difference that now Xk is j-independent as in [141] . To a good approximation, one has 



and 8/c = a^/pxo H — k 2 , (43) 




where the numeric coefficients a,b,b' have the same values as before, given in (]35|) . By substituting j — > t/T, where 
T is the period of oscillations of f(z) ^j, the occupation number of the k-mode becomes 

i ( 2X k VA</> \ . 
n k = CX P I jl * I x interference term . (44) 

As expected rik has a linear exponential growth in terms of t. To compare it with the particle creation in an 
expanding background, we recall that in our case ( a i 2 (c/. (0) and hence (|44|) gives lnrik ~ x 2 for the flat 
background. This confirms the intuition that in general the particle creation via tachyonic resonance is more enhanced 
in a flat background compared to that of an expanding background. This is understandable because the expansion 
of the Universe dilutes the previously produced particles and also reduces the amplitude of the source term 4>(t). 
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One should, however, note that due to the non-trivial dynamics of the interference term, this intuition may not work 
and for some specific regions in the k space expansion of the Universe may enhance the particle production rate. As 
explained in previous subsection, the interference term is time dependent which means the stable or unstable bands 
vary in time. As an example, Figj6] shows cases in which the particle production is more efficient in an expanding 
background. 



III. TACHYONIC RESONANCE IN m 2 <j> 2 /2 THEORY 

The A0 4 /4 inflationary model has the conformal symmetry and the effects of the expanding background can be 
incorporated by the conformal transformations (4>,x) ~ ¥ a ( ? ?)( ( / , j x)- This trick does not apply to m 2 4> 2 /2 model and 
the effect of expansion should be taken care of accordingly. The analysis of the tachyonic resonance for m 2 4> 2 /2 model 
in a flat background has been studied in [l4[ with some brief discussions on the effects of expanding background. Here 
we study the tachyonic resonance for m 2 </> 2 /2 inflationary potential in an expanding background in more details and 
demonstrate that it can have very nontrivial consequences. 

We start with the potential 



V(<f>) 



to , 2 a , 2 A 4 



(45) 



with A > a 2 /2m 2 [14|. The equation for the production of x particles obeys 



Xl + {A k + 2qcos 2z) X k = 



(46) 



where \ — a {t) 3 ^ 2 X ancl 



rat 



71" 

2z ~r 



A, 



4fc 2 



4 and 



2cr$ 
m 2 (J 3/2 



go 

7,3/2 



(47) 



Here prime denotes derivatives with respect to coordinate z and Aq and go indicate the values of the corresponding 
quantities in a flat background. For the quadratic potential, the background Universe during preheating evolves like 
matter domination and a(t) = ao(t/to) 2 ^ 3 and in our conventions we choose do = 1 and to = 1. 

One can incorporate effects of expansion in the tachyonic resonance analysis as follows. The tachyonic regions of 
(|46p are centered around the minimum of w 2 (z), at tj = (j — l/2)n with the scale factor a_(tj) = ((j — l/2)7r) 2 / 3 
whereas the non-tachyonic regions are centered around the maximum of u> 2 (z), at tj = (j — l)ir with the scale factor 
a + (tj) = ((j — l)7r) 2 / 3 . Following the same steps as in section Hi Al thc occupation number is given by 



l^| 2 =exp[ £ 2 *k I! (2cose^). 



(48) 



where is defined as the last non-tachyonic oscillation after which tachyonic resonance starts 



r 1 f A k\ 1 i 



(49) 



where [z] represents the integer part of z. 
Using results of [3] , one finds 



X> = - 



<* A l 



2a- 



<.Al 



1 >2aV qo 



qot 



5/6 



1/2 ' 



(50) 



in which a = 0.85 and tj = n(j — 1/2). Furthermore, the phase accumulation is given by 



9o 



a + {t 3 f 2 



a + b- 



2q a+(tj) 



1/2 



C 1- 



2q a + (t, 



a/2 



In 1 - 



4 



2q a + {t, J 



,1/2 



(51) 



where a = 1.69, b = 2.31 and c = 0.46 [3. 
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FIG. 7: The logarithm of the occupation number as a function of time for A® — 20 and qo — 13. The wavy (blue) curve shows 
the full numerical solution. The (red) dashed curve shows our analytical result (|53p . As expected from (|54p . there is no particle 
production after t = 23. 



Using the approximations in Appendix. [B]for the harmonic sums, we obtain 



E 2X i - v^( (j7r)1/2 - ( ^ )1/2 ) - ( {J7T)1/6 - {j * 7r)1/6 ) 



12a A° k 
Qo 



(52) 



j'.+i 

Finally, by substituting jir — > t after j oscillation, the occupation number as a function of time is 



exp 



— Vmlt 1 

TT \ 



/2 



,1/2 



l2aA l 



.1/6 _ ,1/6 



f[(2 cos 



(53) 



where O k is given by (jSTj) and = j*ir. 

Because of the expansion of the Universe, both q and decrease as t increases, but A^ decreases faster. On the 
other hand, the WKB approximation is valid for 2q — A > 2^/q which quickly reduces to 2q > 2^/q. This condition 
is not satisfied for q < 1. As a result, the expansion of the Universe spoils our approximation and that is why the 
analytical solutions are not as precise as in the case of A</> 4 /4 theory. As one can see from the numerical results in 
Fig. [71 particle production is stopped after some oscillations because of the expansion of the Universe. 

In order to interpret this effect note that by reducing A^ and q, the solutions of Mathieu equation converge to the 
stability bands. Therefore, in the expanding background when the A^ — q curve crosses the stability bands, the particle 
production switches off. From the stability/instability charts of Mathieu equation [h| one finds that preheating ends 
when 



go 

^cnd 



0.8126 1 - 



,4/3 



(54) 



As a results, for initial conditions with A® > 2qg, there can be particle production only if t* > i cn d- This does not 
take into account the effects of back-reaction which will be studied in the next section. 

As already mentioned, expansion of the Universe dilutes the condensate of the produced particles as well as reducing 
(f>(t) as the source of resonance. As a result the resonance in an expanding Universe is expected to be less efficient than 
the non-expanding background. More investigations shows, however, that for some specific modes k the expansion 
of Universe can actually enhance the preheating! There are two effects which can enhance particle production in an 
expanding background. First and less important is the k 2 /a 2 term in preheating equation (|46[) . The expansion of the 
Universe with the effect of k 2 — > k 2 /a(t) 2 can reduce the energy cost of producing particle (but there is a trade-off 
between this effect and the effects of diluting (^-condensate and the reduction of the source term). Second and the 
more important is the effect of varying interference term in an expanding background. 

We first focus on the former effect. The occupation number in a flat background is given by [l4| 



exp 



2j 2a 



(2 cos 0k 



(55) 
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FIG. 8: Left: the logarithm of the occupation number as a function of time for A® = 230 and qo = 125. It shows that the 
expansion of the Universe can enhance the particle creation by reducing the energy cost of particle production via k 2 — » k 2 /a 2 (t). 
The solid curve (blue) shows our full numerical result in an expanding Universe. The lower straight solid line (green) shows 
numerical results for non-expanding universe. The dashed curve (red) shows our analytical results in an expanding background 
(|53[) . Right: logarithm of the occupation number for A\ = 260 and go = 125. This figure shows that because of the time 
varying nature of the interference term, the solutions can escape from the stability bands. The solid curve (blue) represents our 
full numerical result in an expanding Universe. The bottom solid line (green) shows the numerical results for non-expanding 
Universe. The dashed curve (red) shows our analytical solution in an expanding background (|53[1 . 



where 



2q \ 2q) \ 2q 



(56) 



As one can expect from the Floquet theorem, the occupation number has a linear exponential growth with time (or 
number of oscillations) in a flat background. Compare this with our result (|53[) where the occupation number grows 
with an exponent which scales like t 1 / 2 or j 1 ^ 2 . In general this leads to the conclusion that at large t the tachyonic 
resonance is less efficient in the expanding background as compared to the flat background. For the intermediate 
times the situation could be different. To see this note that the second term in the big bracket in (|53|). which 
suppresses the particle creation, scales like t 1 ^ in an expanding background whereas it scales linearly with t in a 
flat background. This difference in scaling results in an enhancement of particle creation for some certain modes in 
an expanding background. However, after some oscillations the second terms in brackets in (|53j) and (|55p become 
negligible compared to the first terms in the corresponding brackets. Approximately this occurs for t ~ 30i* and 
particle production in an expanding background becomes more and more inefficient compared to the static Universe. 
One can see this feature for qo = 125 and A® = 230 in Figj8[ showing that the particle production in an expanding 
background is more efficient for about first 60 oscillations. 

As mentioned above the interference term can play an important role in relative enhancement of particle production 
in an expanding background. In the stability /instability chart of the Mathieu equation there are vast regions where 
there is no particle production due to destructive interference term. However, as we described at the end of subsection 
III BL in an expanding background the phase term is not static and is varying with time. As such and as time advances, 
the solutions in the stability bands can escape from the stability bands. Inside the tachyonic region, far from Ak = 2q 
line, the effect of phase term is suppressed, but by going towards the line Ak = 2q the stable regions are formed. So, 
in a flat background there is no particle creation in this region but as described above, in an expanding background, 
the solutions can escape these stability bands, leading to an enhancement in particle creation (see the right graph in 
Fig®. 



IV. BACK-REACTION OF x-PARTICLES AND END OF PREHEATING 

So far we have neglected the back-reactions of the produced particles on classical background <f> and on \ particles 
occupation number. There are important effects which can change this simple picture. We classify these effects by 
the form of their interactions. 

First, trilinear interaction at one- loop Hartree approximation contributes as a source term in equation of motion. 
This leads to a non-zero vacuum expectation value for 4> which is expected, because as one can see from (JIJ, the 
minimum of the potential is located at a non-zero value of 0. Second, the self-interaction of ^-particles can increase 
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the effective mass of x-particles. This effect can make x _ P ar ticles so heavy that they can not be produced through 
interaction with (f>, terminating the preheating. 



A. Back-reaction in A</> 4 theory 

Let us first ignore the effects of decay of inflaton field through trilinear interaction and consider the inflaton field 
as a background field. The preheating is complete if all the energy from the background field (j) is transferred into 
created ^-particles. One can estimate the energy in the \ particles as 

The 1/a 2 factor in front of the integral has appeared because we are working in the conformal frame and that 
Pconf = Pcomoving/a 2 - Note also that n£ = n^/a 2 (a;), with n£ given in (|36p . The integral over k is cut off at 
kmax ~ ^ L P 2 P x - This upper bound comes from the fact that for a given time x particle creation starts for x > x*, 
which in turn, recalling (|28[) . implies the bound mentioned above. 

In the conformal frame, at one loop Hartree approximation, the dispersion relation is 

^k = V (^ 2 + P xf(x) + 3j^Pj , (58) 

and for our estimate of p x at this stage, we drop the last term and justify this approximation later on in this 
subsection. Moreover, for the period when the particle creation is more pronounced one may also drop the k 2 term 
and approximate w 2 by \(p 2 pxf(x) or 



With this approximation 



where 



uj ~ y/ aa{x)\(p{x)\ . (59) 



1 f d 3 k 1 



is the total number of \ particles produced. In order to perform the above integral we note that it could be written 
in the form 

n x « e -^-^{\v 2 Y' 2 J (Ikk 2 c , (62) 

where C = 4a/(3T) ~ 0.5, A = (126' - 4a)/3pT and B = Ab'^/x/^/pT. This integral cannot be calculated analytically 
and one should approximate it. Comparing the two terms in exponent of the integrand 

Ak 3 (3b' - a) k k .. 

-5-2 = oy 7= ~ °- 7 ^ > 63 



we see that for the range of integration k < y/px one can neglect term Ak 3 and just take —Bk 2 term, reducing the 
integral to an incomplete Gaussian integral. For large values of x, which we are interested in, i.e. for ^/px 3 ' 2 > T/b' w 
2.6, the integral can be computed leading to 

n x « i (X^) 3 / 2 x- 3 '^^ /2 . (64) 



3 Here a = 2.72, as given in (135 b . and should not be mistaken with the scale factor a(x). 
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Preheating completes at a; cop , where 



Px(£co P )a 3 (2; C op) ~ p\ = ^ 4 . (65) 



which happens when 

C ° P Ap5/4 ^ T J V 3 M P Ap5/4 ■ 

As one expects by increasing p preheating shuts off sooner. This is reasonable, since the bigger the value of p, the 
stronger is the trilinear interaction which results in a more efficient x particle production and a stronger back-reaction 
effect. As an example, with p = 0.05 and A = 10~ 14 one finds x CO p — 52.5 which is about 7 oscillations. 

As the density of x particle grows and is seen from (|58p , the effective mass of x particles also grows and production 
of them becomes more costly and eventually terminate the preheating. This happens at x±, and at first loop Hartree 
approximation it is when 

3X'(x 2 )\ Xi ^Xp Xl f. (67) 

On the other hand, the produced ^-particles also back-react on the dynamics of the inflation field <p(t) and may cause 
preheating to stop before completion, making our preheating model inefficient. To check when this can happen we 
note that back reaction of x _ P ar ticles at one loop Hartree approximation level, modifies the (p equation {5} to 

^ + A^ + ^<* 2 >=0. (68) 

The back reaction of x-particles on if becomes important at x 2 , when 

<j(a(x)( X 2 ))\ X2 ^2Xf> 3 . (69) 

Whichever of the two stoping mechanisms happens first marks the end of preheating. We will denote this time by 
a^cop, where x cop — min(x\, x 2 ). The condition of having a successful preheating is then x cop > x cop . 
To evaluate the x\ and X2 we note that within our approximations 

[X } J (2tt)3K| | W | ' 



with \ui\ given in (|59[) . This leads to 



" x(Xl) ~ 3A 7 (^ Ax i) 3/ V > ( 71a ) 
n x (x 2 ) » 2A (p A x 2 y 1/2 <fi 3 . (71b) 



Using n x given in (|6^|) we end up with 



x -m e ^r „ " ^ 3/2 pVi ^ 2|5i p3/4; (72a) 

(72b) 



A \T J y ~ Ap 5 / 4 ' 

As we see x\ depends on A' whereas x 2 depends on A. For A ~ A' ~ 10 _ and with p = 0.05, we find that 
x\ ~ x 2 — 52.5. This means that for this choice of parameters ir cop ~ 52.5 which is very close to the completion of 
preheating time x cop and hence we expect an efficient preheating. 

Given the above expressions one may ask for which range of parameters p, A, A' the "efficient preheating" condition 
x cop > x cop is satisfied. The comparison between x 2 and x cop is somewhat straightforward, noting that in the range 
of parameters that we are mainly interested in, n x /uj is a monotonic function of x and hence if a(x 2 ) < 8 (or 
equivalently x 2 < 55) x 2 < x cop and vice-versa. Moreover, one can argue that either x\> U > x\ with X\X 2 > U, 
or x\ < U < x\ with 3:1X2 < U, where U = 6A'/(Ap 2 ). In the former case, the x cop > x cop condition roughly boils 
down to 500p 2 jr > 1. In the latter case, when x cop = Xx, one may again show that a similar bound holds. This 
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condition together with ([3]) specifies the range of parameters for which we have a simple slow-roll A(/> 4 inflation as well 
as efficient preheating: 

2 x 1CT 3 <p 2 ^ < 1 . (73) 
A 

For typical A ~ A' values, this leads to 2 x 1CP 3 < p 2 < 1. 

It is also notable that aa(x)(x 2 ) ~ V°" a (f)_" x. term shifts the minimum of the potential of <p from w — and we 

expect the end point of preheating to be not far from the <p m i n . As a rough estimate of this minimum value (assuming 
that x\ ~ x 2 ), at this approximation level, is <p 2 min ~ a 2 (x2)cr 2 /(6AA'). For our estimates, where a(x 2 ) ~ 8, this rough 
estimate is not far from the the global minimum of the potential (fTJ) which is at 4>o = cr 2 /(2AA'). 

Before ending this subsection let us briefly discuss the production of cp particles and re-scattering of x-particles 
in the theory with potential (Q]). At tree level the (£k modes will be sourced by a term like X(ip 2 }. This term being 
positive and periodic, following the discussions of [l7j, leads to stochastic resonant production of </?k modes, which 
at large x, has an average exponential growth in x. This effect compared to the tachyonic resonant production of \ 
particles is very small and one may safely conclude that during preheating mainly x-particles are produced. 

At one loop level there is a contribution to the equation of motion for the %k mode proportional to \a 2 <p 2 /k 2 oc 
A 2 p 2 (-^p) 2 f 2 (x). This, unlike the tachyonic source term, pxf(x), is always positive. Nonetheless being proportional 
to X 2 p 2 , this one loop contribution, is too small compared to the tree level pxf(x) term. With the trilinear 4>x 2 
coupling, produced x particles will not back-react on the the production of ip^ modes, at first loop level beyond the 
Hartree approximation, which we have already discussed. This is in contrast with the four-leg 4> 2 x 2 interaction, where 
a non-zero (x 2 ) contribute to the equation of motion of ip^. 



B. Back-reaction in m -theory with trilinear interaction 

The important point in this case is that following and discussions leading to it, we note that in the m 2 2 
theory there is a definite time, i cn d, after which there is no particle production at all. Therefore, to have an efficient 
preheating the energy transfer form the <fi background into x particles should happen before this time. 

As in A0 4 /4 case lets us estimate the time at which the energy transferred into the x particles becomes comparable 
to the background energy. The energy density of x-particles after time t is 



Px 



d 3 fc 
(2^)3 



KK = 



1 



d 3 fc 



a 3 {t) J (2tt) 



\u) k \rij , 



where rij is given in (|53p and from potential (|45[) one reads that at one loop Hartree approximation 



a<p(t) 



3A 



a?{t) (a(t)) 3 / 2 a 3 (t) 



+ ^7V(X 2 > = k 2 r^ + «rt-V(«) + 3At~ 2 (x 2 ) 



(74) 



(75) 



with tp(t) = $0 sin mi. The main contribution to the integral (|74p comes from the period when the tachyonic resonance 
is at work. In this period one may drop k 2 /a 2 term in (J73J). For the current estimation we also drop the back-reaction 
of x-partideSj to which we will return later. Therefore, 



Inserting (fT6) into (f?4")) we find 



where 



Px 



d 3 fc 



(76) 



(77) 



d 3 fc 
(2tt) 3 



exp 



7T \ 



- t 



1/2 



12aA° 



(78) 



In the above, as in the A</) 4 case, the upper bound on k, k max , comes from the fact that particle creation becomes 
efficient for t > i* = nj* with given in (14^1) . that is k max < m^f^-t x l^ '. To perform the integral we note that it has 
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the same form as in (|62|) . with a cubic and a quadratic term in the exponent, where now A = f^^ 1 , B = ^f^- t 1 / 6 

and hence Ak/B < 2/3. One may again drop the Ak 3 term and approximate the integral by an incomplete Gaussian 
integral: 



dfcfc 2 e- Bfc ~ ? ° ; J -'/ 4 e^ (79) 



2vr 2 y 512(3a) 3 / 2 



2 



where £ = 8a /ir ~ 2 and in the last step we have approximated the integral for qo^ ^ 1- 

Preheating is complete at t cop when the energy of x-particles becomes comparable to the energy in the the inflaton 
condensate at the beginning of preheating in the same comoving volume, i.e. 

Px a\t) =(n x H) -Jm 2 * 2 . (80) 

Using the expressions given above we obtain 

« 256V2(3a) 3 / 2 9o " 5/4 ~ 1474 q~ 5/4 

For qo = 50 and <I>o/m ~ 10 5 , the above is satisfied for t cop ~ 3.5 which is slightly more than one period time (which is 
t = tt). In order to have an efficient preheating, we should demand that the preheating completes before the particle 
creation stops. The latter happens at i e nd given in (|54[) . i cn d > t cop imposes a lower bound on qg. For &o/m ~ 10 5 , 
we find 

qo > q c - 13 . (82) 

Noting that qo = 2<j$>o/m 2 , this can be used to impose a lower bound on a, the scale involved in trilinear coupling, 
cr/m > 6.5 m/$ — 6.5 x 10 -5 . 

As (x 2 ) increases particle effective mass increases, cf. (|75l) . which in turn can stop preheating. One needs to also 
verify that this happens after the completion of preheating, that is 

3Ai~ 2 (x 2 ) < <rt _ V(*) (83) 

at t cop . Noting that 

(X 2 )^n x ±- (84) 

\LU\ 

with n x given in (79]), X(x 2 ) ~ at(p(t) happens when (q ty 3/4 e c ^ ~ -|x(^) 3/2 . The condition (|53J) is satisfied 

if A < 0.17(^) 2 q 2 , = 0-67 (^) 2 . Noting the condition for positivity of the potential 1)45 p. an efficient preheating 
scenario happens if A is in the very tight range: 



O.sf-") 2 < A<0.67f- N ) 2 . (85) 



For o-/m - 10~ 4 (when g = 20) that is, 5 x 10~ 9 < A < 6.7 x 10" 9 . 

In the m 2 4> case, with the potential (|45[) . at one loop Hartree approximation level the equation of motion for the 
background cp field is modified as 

<P + ™ 2 <p + ^7^(X 2 ) = , (86) 

where a 3 / 2 = t, <p = a 3//2 $ and x = a 3 / 2 %. Using (|84|) and (|76|) . we observe that the force term in (|86|) at the 
completion of preheating time can vanish for a non-zero (p m in- rn 2 (p m i n ~ cri _1 n x /2ii;| 4 _ 4 which upon using (|80j) 
we obtain 

m v? mi „ ~ at * = (87) 

or fmin ~ 5^0- It is interesting that our rough estimates and considering the one loop Hartree approximation 
reproduces the order of magnitudes of the minimum obtained in numerical analysis of [l4| (see FIG3 in Ref . [lij ) . 
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V. SUMMARY 



In this paper we have studied tachyonic resonance preheating via trilinear interaction in an expanding background. 
Because of the three-leg interaction the conformal symmetry in A0 4 /4 inflationary theory is broken. This induces 
a non-periodic source term in resonant x field equation. Interestingly, one observes that the particle creation has a 
non-linear exponential enhancement with the leading exponent ~ x 3 ^ 2 . This is in contrast to particle creation via 
parametric resonance from the four-legs interaction in this theory which preserves the conformal invariance and a 
linear exponential growth of particle creation is obtained. Besides the non-linear exponential growth, the interference 
term obtained from the accumulation of the phase term in non-tachyonic scattering regions has a very non-trivial 
behavior in an expanding background. It is shown that there are "semi stability bands" where djQ 3 ^ — cosO^ = for 
some oscillations j s . As a result, the particle creation for the corresponding modes are highly suppressed. However, 
due to time varying nature of the phase term, the semi-stability bands are washed out as time goes by. 

The tachyonic resonance preheating in an expanding background for m 2 (j> 2 /2 theory was also studied. In general 
the expansion of the Universe suppresses the particle creation. However, due to time varying nature of the interference 
term, the expansion of the Universe can actually enhance the particle creation for certain modes. 

We studied in some detail, the back reaction of preheat x-particles on the dynamics of preheating and the background 
inflaton field <f>. As we argued demanding an efficient preheating (i.e. demanding that particle production ends not 
before most of the energy in the background ip field is transferred into the x-particles) imposes strong bounds on the 
parameters of the potential in both A0 4 and m 2 (f> 2 cases. 

For typical values of the parameters of the potential for X(j) 4 theory preheating is complete and lasts for < 10 
oscillations while for m 2 (j) 2 in few oscillations the energy transfer to x _ P ar ticles seems to be complete. 
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APPENDIX A: PROPERTIES OF JACOBI ELLIPTIC COSINE FUNCTIONS cn(x\m) 

Here we briefly review the properties of the Jacobi elliptic cosine function cn(z\m). Note that we follow the 
convention of 1101] which differs from the convention of [17| . A useful representation of the Jacobi elliptic cosine 
functions is [18| 

/ , , 2tt ^ q(m) n+1 / 2 f nz \ 

cn(z\m) = 1 =— , > / , , cos (2n + \) n „, r , (Al) 

v 1 ' y/mK{m) £rj,l + q{mf n + l \ K ' 2K(m) J y ' 

where 

nK (1 — to 



q(m) = exp 



K{m) 



and K(m) is the complete elliptic integral of the first kind. 

The Jacobi elliptic cosine is periodic with period T = 4 AT (to) and has zeros at 



xj = (J \)T, xj = (j + \)T. (A2) 



The minima and maxima are at 



xj ax = (j - l)T, *7 n = (J - ~)T. (A3) 



For the special case of m = 1/2 this formula reduces to 



1 2V2^ e-^+Vs) 



f(x) = cn(x\=-) = ^Ly — ^-rvcos (2n+l)—4jT } (A4) 
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As can be seen from above, f(x) is even and periodic with the periodicity T = 4K"(±) = r "^/ 4) « 7.416. 

APPENDIX B: APPROXIMATIONS FOR GENERALIZED HARMONIC NUMBERS 

Here we present some useful formula for the sums of the harmonic numbers which are used to calculate the sum 
over j for the occupation number in (|25| , (|33| and (|52|) . 
With straightforward algebra one can check that 

£^4 s/2 +?" 2 4 < bi > 



±1, 2j ^ + y^i m 



i=l 

and 



£V2^I^^V 2 (B3) 



which arc numerically confirmed to have a very good accuracy. 

One can generalized these approximations for arbitrary power s ^ -1 



y i< ~ — + 1/ + . (B4) 



APPENDIX C: APPROXIMATION FOR CALCULATING f x i ax _ T/4 y/xf{x)dx 



In this section we would like to examine the approximation used in calculating the integral in ([23]): 

xf a* +T/i + T 

y/xf(x) dx ~ yjx?** / t y/f{x)dx. (CI) 



-T/4 



Defining a; = £™ a2: + u, one finds 



f^+T/A ,T/4 _ 

y^TRdx = Jxf^ / /i + -^Jf(xT ax + u ) d * 

xf ax -T/i v ./-T/4 V j 



T/4 



^ / ./i + ^V/Rd"- (C2) 



J / , \ max 
-T/4 V ^ 



In the last step we have used the relation /(xj 111 * + u) = f(x] (see Appendix |A")) . Expanding the first square root in 
the integral and noting that f(x) is an even function one finds 



For our case, we need j > 1, so the first correction can lead to maximum 6% error. Since the higher corrections are 
inversely related to square of j — 1, the errors decay very quickly. For example for j = 4 the error is less than 1%. 
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APPENDIX D: APPROXIMATION FOR 



CALCULATING fl W Jr + f(x)dx FOR < r < 1 



Here we demonstrate the approximations used to calculate the integral in (f3"Tj) and (|3"2"|) 
From [ATI one finds that 

cn(x + xj' ax \m) — cn(x\m) and cn(x + x™ m \m) = — cn(x|m) (Dl) 

The function r + / (a;) with < r < 1 has zeros at 



a;7(r) = xf ax {r) + K arccos(-r)|- = xf m (r) - 2K - )+K[ arccos(-r)|- 



4(r) = a^f (r) - K (arccos(-r)|0 = xf n (r) + 2K 



- X ^arccos(-r)|i ) (1)2) 



where K{x\m) is the incomplete elliptic integral of the first kind with parameter m and amplitude x We propose 
following approximations: 

x~ (r) 



_r m 

with a = ,/f — W- = 2.72 , b = 3.748 and 



Vr + f(x)d x~a + br, (D3) 



x+(r) 



y/\r + f{x)\dx ~a-b'r, (D4) 

with V = 2.864. The fact that b 1 > has the important consequence that the occupation number ^ is more 
suppressed for higher values of k, see subsection III Bl for the details. 

Both of these approximations are verified numerically to a good accuracy. Of course, to have higher accuracies one 
can keep higher powers of r in above expansions. 
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